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We present a theory of a single-electron transistor exchange-coupled to a localized spin. We show 
how to gain detailed quantitative knowledge about the attached spin such as spin size, exchange cou- 
pling strength, Lande g-factor, and spin decay time Ti by utilizing a robust blockade phenomenon 
of DC magnetotransport with accompanying noise enhancement. Our studies are of particular rel- 
evance to spin-resolved scanning single-electron transistor microscopy, electronic transport through 
nanomagnets, and the effect of hyperfine interaction on transport electrons by surrounding nuclear 
spins. 

PACS numbers: 72.25.Rb, 72.70.+m, 73.23.Hk, 73.63.Kv 



The selective detection and manipulation of single 
spins are of fundamental importance for the realization of 
qubits as basic building-blocks in quantum computers 1 . 
In a solid-state environment, spins experience relaxation 
and decoherence which limit the functionality of quan- 
tum information processing. Therefore it is necessary 
to shed light onto the time scales of those effects. In 
addition, properties such as the coupling between spins 
as well as to an external magnetic field are important. 
Recently, the electrical time-resolved read-out of an indi- 
vidual electron spin enabled the measurement of the re- 
laxation time^ T\. Furthermore, Wabnig el al. have pro- 
posed obtaining the T\ and T 2 times by high-frequency 
(GHz) noise measurements^. Such experiments would 
be very challenging and we suggest here an alternative 
route to access at least T\ which avoids time-resolved or 
high-frequency setups. 

In this letter we explore the nonlinear current and 
low-frequency noise of a single-electron transistor (SET) 
where the spin of the quantum dot (QD) electron is cou- 
pled via the exchange interaction to a second localized 
spin. We find that tunneling spectroscopy in a constant 
magnetic field reveals a robust current blockade region 
due to population trapping that allows for the determi- 
nation of T\ in a simple manner. The proposed setup 
could be used as a spin-resolved version of the scan- 
ning single-electron transistor (SSET) 4 , as alternative 
to single-spin detection by magnetic resonance force mi- 
croscopy (MRFM)£. 

There are two further situations for which our work is 
of relevance: Electronic transport through nanomagnets 
(e.g. QDs doped with Mn-ions^) with an attached spin of 
5/2 and an anisotropic exchange coupling; and as a start- 
ing point for studying the effect of hyperfine interaction 
on transport electrons by surrounding nuclear spins in 
QDs (e.g. as shown in Refs. 7,? for double-QDs). 

Our model consists of a SET with constant lead 
couplings TL t R and a QD electron spin S 2 which is 
isotropically exchange-coupled to an additional spin 
Si with strength J as sketched in Fig. QJl. The corre- 
sponding Hamiltonian (for the exchange term see also 
Refi^), essentially an Anderson Model with exchange- 



coupled spin Si, reads H 
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Hqd — (£ini + mB-S i )+JS 1 -g 2 , (1) 

i=l,2 
ka<72 

with £j the single-particle energy of the QD (i = 2) 
and the attached spin level (i = 1), rii = S a .d ai - 
Throughout this work the SET operates in the Coulomb 
blockade regime, i.e., only a single electron can enter 
the SET. This restricts the dimension of the accessible 
Hilbert space. The attached state with spin S\ is al- 
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FIG. 1: (color online), (a) Sketch of the single-electron tran- 
sistor (SET) with lead couplings Vi, : r and electronic level 
with spin S2 and attached spin Si via exchange coupling J. 
(b) Eigenenergies of -Hqd (singlet S, triplets T,T±, unoccu- 
pied SET-level V±) vs. magnetic field B for (71,2 > (see 
Tab. U]). Arrows indicate allowed transitions due to single- 
electron tunneling, (c) Steady-state SET current vs. sym- 
metric bias voltage V and magnetic field B; Solid lines: cur- 
rent steps. Dotted lines: blocked transitions; Shaded region I: 
current blockade, bunched electron transfer (see also Fig.[2ji). 
Region II: see text. 
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TABLE I: Eigenenergies and -states of Hqy> (0 for S = Be z , 
el 1 ' = figjB/J, e = £i + £ 2 , £ z s e! 1 ' + £i 2) , Ae 2 = £i x) - 
£i 2) , d/ 2 s yl + c 2 /4 , c 3/4 = Ae z ± c, c = V4Ae| + 1. 

(All energies are scaled by the exchange coupling strength J 
throughout the paper. Realistic values for J can be estimated, 
e.g., by the dipole-dipole interaction between two electron 
spins^ J « (fo/f) 3 0.1/xeV with distance r and ro = lnm. ) 

ways occupied with one electron. The Zeeman-energies 
ptgiB • Si contain the Lande factors gi and magnetic field 
B. J denotes the exchange coupling strength. The op- 
erators dj/di (i = 1,2) are the creation/annihilation op- 
erators in the attached spin/SET, and ct aa ./cf~ a(T are the 
creation/annihilation operators in lead a = L, R for mo- 
mentum k and spin a. 

For the sake of simplicity and clarity we consider a 
spin Si = 1/2 throughout this work. The corresponding 
eigenenergies and -states of Hqv are presented in Tab.Q] 
with the magnetic field-dependent spectrum shown in 
Fig. QJ>. Note that the coefficients Cj (z = 1 ... 4) depend 
on the magnetic field. The Coulomb interaction between 
the QD electron and the side-electron simply produces 
an overall-shift of the spectrum and is therefore not con- 
sidered further. 

We treat the coupling to the contacts in second- 
order perturbation theory (sequential tunneling ap- 
proximation). The standard Born-Markov-Secular 
approximatio n ' for B > (non-degenerate spec- 
trum) leads to rate equations in the energy-eigenbasis. 
They can be cast into the compact form p — Cp with 

P = (P(T0) > P(io) > Pr + j PT- j PTi Ps) and C = C L + C R . 
The Llouvillians (6 x 6-matrices) Clir contain the cou- 
plings T a L/R = 2iTj2 k \t k L/Ra\ 2 S(s - e kL/Ra .), the Fermi 
functions /^(ei) for the occupation of the L/i?-contacts 
with the excitation energy ej and respective bias volt- 
age ±V/2, and the Clebsch-Gordan coefficients c^/ci, 
1/ci, C4/C2, I/C2 (see Tab. |T|. Counting fields e ±lx are 
included in the right lead Liouvillian Cr = £r(x) 13 - 
The cumulant generating function is given by S(x, t) = 
In |Tr[exp(£(x)<)p]} with steady-state density matrix p 
given by C(0)p — 0. The cumulants are obtained via 
Cfc(t) = (-i) k d^S(x,t)\^ x=0 y The steady-state current 

is (/) = eC\ and the zero-frequency Fano factor (DC 
noise) is defined as F = CijC\. The Fano factor in- 
dicates sub(super)-Poissonian electron transfer if it is 
smaller (larger) than unity. 
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FIG. 2: (color online), (a) Fano factor vs. symmetric bias 
voltage V and magnetic field B, parameters as in Fig. []J; 
region P: Poissonian charge transfer; Inset: configuration 
space with single-electron tunneling transitions for region / 
in Fig.[TJ;, Transition \V-) — > \S) is exponentially suppressed, 
(b) Fano factor vs. spin decay rate 71 for symmetric coupling 
r = \V L / R and various coefficients C4 oc B in region / (see 
Fig. [2^,); Inset: Configuration space with indicated spin flip 
rates. 



The spectrum of -ffQD (Tab.|TJ itself is not directly ac- 
cessible by transport spectroscopy. Rather, by applying 
a bias voltage V to the SET, the excitation spectrum is 
probed. The six transitions with nonvanishing Clebsch- 
Gordan coefficients are indicated in Fig. [TJa. Fig. Q}j 
shows that the current vs. magnetic field B and bias volt- 
age V is a direct map of those transitions. The six steps 
correspond to excitations in the eigenspectrum (assuming 
the temperature is low, otherwise the transitions would 
be smeared out). Between steps the current is nearly 
constant for fixed, but changes with variable magnetic 
field due to the dependence of the Clebsch-Gordan coef- 
ficients. At zero-magnetic field, two current steps occur 
corresponding to tunneling transitions from the unoccu- 
pied states to the singlet- and degenerate triplet states. 
The bias difference of the steps is 2J. From the posi- 
tions of the current steps one easily obtains the g-factors 
of both spins by fitting the excitation energies calculated 
by the differences of eigenenergies in Tab.|H We point out 
that the specific dependence of the excitation spectrum 
on the magnetic field can appear very different when one 
of the g-factors is negative (as widely occurs in confined 
electron systems) e.g. for gi < negative differential 
conductance appears for certain transitions. 

Remarkably, certain transitions are completely miss- 
ing (dotted lines in Fig. [TJ;) even though the correspond- 
ing Clebsch-Gordan coefficient is not zero. This is to 
be contrasted with spin blockade for which the relevant 
coefficients disappear 5 ^. This blockade here is caused 
by the topology of the configuration space. Let us fo- 
cus on the lowest lying transition in Fig. [TJ; where the 
configuration space consists of three states: \V+), \V—), 
and I S) in the inset of Fig. [2^. Since the transition 
energy for \V—) — > \S) is too large to be supplied by 
the lead electrons in that bias range, the rate for that 
process is exponentially suppressed. Consequently, the 
system gets trapped in the state \V-) once it reaches 
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there, the spin Si is completely polarized and the cur- 
rent is blocked. Such population trapping effect can be 
expected in any system with a configuration space with 
one or more singly-connected states. We emphasize that 
the current blockade persists for an arbitrary choice of J, 
gi,2, Si, and even for ferromagnetic leads T^ L R ^ T^ L R . 
However, the blockade will be removed with increasing 
temperature T since the transition |V_) — > |S) becomes 
more likely (compare inset of Fig. [2^). In regions I 
and II the dynamics can be described by an effective 
three-state model assuming low temperatures such that 
/ii( e si) — Jr{^s2) — and excitations involving triplet 
states are negligible: 



P\10) = {-^Lh(e s i)p\io) + Fl/l ( £ si) + T R ]p ls) }c 2 , 

c 2 

Ptfo) = {-ri/i(e s2 )p| T0 > + Pi/i (e s2 ) +r R ]p| S> }-|, 

P\s) 



same magnetic field (assuming 71 <C T so that (/)// does 
not depend on 71): 



{F L fL(e s i)p\io) + c 2 4 T L f L (e s2 )p U0) 
[Tl\JE + +TbcI]p\s)W (2) 



with e sl/s2 =e 2 ±ei 1) /2-l/4-c/2and/^(-) = l-f L (-). 

The Fano factor in region I is Fi = l+2c!r fl /(r L +r_R) 
and clearly indicates super-Poissonian electron transfer 
since C4 > 1 for B > (lower left corner of Fig. [2^) . A 
very similar expression has been derived in Rcf. 15 with 
the interaction parameter <I> replacing C4. The underly- 
ing mechanism is related to thermally-activated bunching 
of tunneling events^: as long as the system is in state 
|V_) no tunneling occurs. However, with an exponen- 
tially small probability, thermally-excited lead electrons 
can enter the state |S) and a small bunch of tunneling 
transitions |V+) |S) may take place. Outside of region 
I the noise is Poissonian (P) or sub-Poissonian as can be 
seen in Fig. [2^l for symmetric lead couplings Tl = T R . 
For asymmetric couplings, it is also possible to observe 
F > 1 in other regions, e.g., in region II when the con- 
dition T L /T R > (c 4 - 1)7 2c 4 is fulfilled. 

We can exploit the current blockade mechanism in or- 
der to measure the spin relaxation rate 71 = 1 jT\ of Si . 
In the inset of Fig. [2Jd we show that the formerly singly- 
connected state \VJ) is additionally linked with state 
|V+) via spin flip transitions with rates 7+ = 71 /(e£ ), 
7_ = 71 [l — /(ei )] • Such a mono-exponential spin de- 
cay can occur, e.g., due to spin-orbit coupling 17 . The 
steady-state current is then directly related to the spin- 
relaxation rate 1/T\ and it can be simply measured by 
the ratio of plateau currents in regions I and II at the 
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with Ai = 4/{TT L 4), A 2 = [l- f(eP)](T L + 2T R ), 

A 3 = (cjr + T L [1 - /(ei 1 ')] + r R ), and r = T L + T R . 

The Fano factor vs. 71 in Fig. shows that by in- 
creasing the spin relaxation the charge transfer turns 
from super- to sub-Poissonian. For large magnetic fields 
the Fano factor approaches unity: for small 71 from 
the super-Poissonian side and for large 71 from the sub- 
Poissonian side. Hence, the zero-frequency noise provides 
a sensitive quantitive indicator of the spin decay. 

In region / the exponential suppression of the 
current can also be lifted by algebraic cotunneling 
contributions 18 . To ensure that the effect of spin relax- 
ation on the SET transport is not obscured by higher- 
order tunneling processes, the rates T^/ R have to be 
small with respect to the temperature. 

The discussed effects are robust against background 
charge fluctuations. We have checked that by consid- 
ering fluctuating energy levels (t) = + 1/ y/n & (t) 
(i = 1,2) with = and (6(t)&(t0) = 5 l3 8{t-t'). 

Second-order perturbation theory in the fluc- 
tuation strength leads to the Lindblad form 11 

£rfP = Si r i~ i n iP n i ~ |{ n < i P}]i which leaves the 
relevant subspace invariant. 

In summary, we have studied a Coulomb-blockaded 
SET with exchange-coupled 1/2-spin in a constant mag- 
netic field. The excitation spectrum (current vs. mag- 
netic field and bias voltage) shows rich structure enabling 
the measurement of the exchange-coupling strength J 
and the Lande g- factors. Due to classical population 
trapping the lowest-lying transition is blocked indepen- 
dent of J, the g-factors and the spin size - the current 
is exponentially suppressed and the electron transfer is 
super-Poissonian. As we demonstrate, this robust phe- 
nomenon allows the direct measurement of the decay rate 
1/Ti of the probed spin by DC current and noise mea- 
surement. Our method can be readily used to implement 
larger Si encountered in nanomagnets or to describe the 
nuclear spin environment of a QD. 
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